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HIGHLIGHTS 


•  We  propose  to  perform  adaptive  SoC  estimation  using  multiple  models. 

•  We  build  an  SoC  estimation  algorithm  via  fusing  estimates  based  on  multiple  models. 

•  The  algorithm  achieves  effective  SoC  estimation  against  complicated  battery  dynamics. 

•  We  validate  the  effectiveness  of  the  algorithm  via  both  simulation  and  experiments. 

•  The  multi-model  framework  can  be  extended  to  other  tasks  for  battery  management. 
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Adaptive  estimation  of  the  state-of-charge  (SoC)  for  batteries  is  increasingly  appealing,  thanks  to  its 
ability  to  accommodate  uncertain  or  time-varying  model  parameters.  We  propose  to  improve  the 
adaptive  SoC  estimation  using  multiple  models  in  this  study,  developing  a  unique  algorithm  called  mm- 
AdaSoC.  Specifically,  two  submodels  in  state-space  form  are  generated  from  a  modified  Nernst  battery 
model.  Both  are  shown  to  be  locally  observable  with  admissible  inputs.  The  iterated  extended  Kalman 
filter  (IEKF)  is  then  applied  to  each  submodel  in  parallel,  estimating  simultaneously  the  SoC  variable  and 
unknown  parameters.  The  SoC  estimates  obtained  from  the  two  separately  implemented  IEKFs  are  fused 
to  yield  the  final  overall  SoC  estimates,  which  tend  to  have  higher  accuracy  than  those  obtained  from  a 
single-model.  Its  effectiveness  is  demonstrated  using  simulation  and  experiments.  The  notion  of  multi¬ 
model  estimation  can  be  extended  promisingly  to  the  development  of  many  other  advanced  battery 
management  and  control  strategies. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Industrial  applications  of  batteries  usually  require  a  well- 
designed  management  system  for  operational  safety  and  perfor¬ 
mance,  which  monitors  the  running  status  and  regulates  the 
charging/discharging  processes  [1].  One  of  its  fundamental  func¬ 
tions  is  to  estimate  the  state-of-charge  (SoC),  i.e.,  the  percentage 
ratio  of  the  present  battery  capacity  over  its  maximum  capacity. 


*  Corresponding  author.  Tel.:  +1  617  621  7500;  fax:  +1  617  621  7550. 
E-mail  address:  yebinwang@ieee.org  (Y.  Wang). 

0378-7753 /$  -  see  front  matter  ©  2013  Elsevier  B.V.  All  rights  reserved. 
http://dx.doi.Org/10.1016/j.jpowsour.2013.12.005 


1.1.  Literature  review 

SoC  estimation  has  remained  an  active  research  field  during  the 
past  years,  and  the  reader  may  refer  to  Ref.  [2]  for  a  survey.  A 
notable  trend  in  this  area  is  the  increasing  emphasis  on  model- 
based  estimation  methods.  The  dynamic  models,  derived  from 
either  equivalent  circuits  or  electrochemical  principles,  facilitate 
the  assimilation  of  the  battery  data  and  lead  to  real-time  SoC 
estimation  with  bounded  errors.  While  battery  modeling  has  been 
well-accomplished  [3],  more  attention  is  being  geared  towards  the 
development  of  estimation  algorithms.  Application  of  the  Kalman 
filtering  (KF)  techniques  has  been  remarkable  in  this  respect.  The 
classical  linear  KF  and  its  extensions  to  nonlinear  systems, 
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including  the  extended  KF  (EKF),  unscented  KF  (UKF),  iterated 
extended  KF  (IEKF),  have  been  used  to  deal  with  SoC  estimation 
based  on  electrochemical  and  equivalent  circuit  models,  see  Refs. 
[4-13].  A  variety  of  other  state  observers  originating  from  control 
theory  have  also  played  a  role  in  constructing  SoC  estimators.  Here, 
we  highlight  the  sliding  mode  observer  [14],  adaptive  model 
reference  observer  [15],  Lyapunov-based  observer  [16]  and  PDE- 
based  observer  [17,18]. 

Since  a  good  model  is  a  prerequisite,  model-based  SoC  esti¬ 
mation  typically  is  after  the  procedures  of  dynamic  modeling  and 
parameter  identification.  However,  accurate  identification  is 
challenging.  First,  the  parameters  in  a  battery  model  are  often 
subject  to  changes  with  time  and  operational  conditions.  For 
instance,  the  internal  resistance  will  rise  and  the  capacity  di¬ 
minishes  as  a  result  of  battery  aging.  Another  example  is  the 
charging  and  discharging  efficiencies,  which  are  dependent  on 
the  SoC,  magnitude  of  current  and  temperature.  Second,  the 
parameters  may  differ  from  one  battery  to  another,  making 
identification  for  each  battery  rather  cumbersome.  Therefore, 
adaptive  approaches  are  more  desirable,  merging  both  identifi¬ 
cation  and  SoC  estimation  in  one  step.  As  shown  in  Fig.  1,  an 
adaptive  SoC  estimator  gives  not  only  the  SoC  estimates  but  also 
the  estimates  of  the  model  parameters  in  real  time  after  assim¬ 
ilating  the  current-voltage  data  on  the  basis  of  a  model.  The 
parameter  estimates  will  then  be  used  to  update  the  model  to  aid 
the  next-step  estimation. 

Adaptive  SoC  estimation  has  been  attracting  considerable 
attention  in  the  recent  literature.  An  adaptive  EKF-based  SoC  esti¬ 
mator  is  designed  in  Ref.  [9],  which  interacts  with  a  parameter 
estimator.  In  Ref.  [11],  state  augmentation  is  conducted  to  incor¬ 
porate  the  SoC  variable  and  model  parameters,  and  then  the  UKF  is 
applied  to  estimate  the  augmented  state.  However,  the  conver¬ 
gence,  and  as  a  result,  the  accuracy,  are  noted  to  be  difficult  to 
guarantee.  In  Ref.  [13],  an  adaptive  SoC  estimator  is  developed  us¬ 
ing  the  IEKF,  guided  by  an  analysis  of  the  observability/identifi- 
ability.  Novel  adaptive  PDE  observers  for  SoC  estimation  have  also 
been  reported  in  Ref.  [19].  It  should  be  noted  that  all  these  existing 
approaches  are  based  on  a  single  battery  model,  and  here  we 
instead  propose  to  exploit  multiple  models  for  better  estimation 
performance. 

1.2.  Statement  of  contributions 

In  this  paper,  we  aim  to  achieve  adaptive,  high-fidelity  and 
easy -to -implement  SoC  estimation.  For  this  purpose,  we  seam¬ 
lessly  link  the  notion  of  ‘multiple  models’  and  adaptive  SoC 
estimation.  A  multitude  of  models,  compared  to  a  single  one,  can 
give  a  better  description  of  complicated  uncertain  dynamics  [20- 
22],  thus  particularly  suitable  to  deal  with  the  tasks  relevant  to 
batteries.  The  design  of  the  adaptive  SoC  estimator  partially 
builds  on  our  previous  work  [13].  In  that  work,  we  propose  an 
adaptive  approach  for  SoC  estimation  via  IEKF-based  simulta¬ 
neous  state  and  parameter  estimation.  While  credible  estimation 
is  observed,  the  accuracy  is  still  limited  in  Ref.  [13]  by  the 


mismatch  between  the  model  and  the  true  system.  This  fact 
motivates  the  development  of  the  MM-AdaSoC  algorithm  in  this 
paper. 

An  overview  of  the  construction  of  MM-AdaSoC  is  as  follows. 
First,  multiple  submodels  are  brought  up  from  a  modified  Nernst 
battery  model  by  fixing  some  parameters  and  assuming  the  others 
unknown.  Each  submodel  is  shown  locally  observable  with  ad¬ 
missible  inputs  by  rigorous  analysis.  Then,  an  adaptive  SoC  esti¬ 
mation  scheme  will  be  implemented  simultaneously  but 
separately  to  each  submodel,  with  the  submodel  in  each  imple¬ 
mentation  assumed  true.  The  SoC  estimates  resulting  from 
different  submodels  will  be  fused  in  the  light  of  a  certain  strategy 
to  obtain  the  final  estimate.  As  such,  we  boost  the  accuracy  of  SoC 
estimation  despite  the  presence  of  uncertainties  plaguing  battery 
models. 

The  main  contributions  of  this  paper  lie  in  two  aspects.  First,  this 
is  the  first  known  study  of  multi-model  adaptive  SoC  estimation  to 
the  best  of  our  knowledge,  and  it  is  shown  that  the  proposed  mm- 
AdaSoC  algorithm  provides  more  accurate  estimation  while 
maintaining  a  good  balance  over  the  computational  cost.  Second, 
we  introduce  the  multi-model  framework  for  battery  management 
and  control.  In  addition  to  the  proposed  MM-AdaSoC,  we  discuss 
various  other  ways  to  enhance  SoC  estimation  using  multiple 
models.  Many  existing  battery  management  strategies  can  also  be 
improved  on  a  multi-model  basis. 

1.3.  Organization 

The  rest  of  the  paper  is  organized  follows.  Section  2  presents  a 
basic  review  of  the  multi-model  estimation  theory.  Section  3  de¬ 
scribes  the  model  construction  and  gives  observability  analysis. 
Section  4  combines  adaptive  SoC  estimation  and  multi-model 
estimation  to  establish  the  MM-AdaSoC  algorithm,  the  effective¬ 
ness  of  which  is  validated  in  Section  5  by  simulation  and  experi¬ 
mental  results.  Finally,  Section  6  gathers  our  conclusions  and  ideas 
for  future  work. 

2.  Basics  of  multi-model  estimation 

The  structure  of  a  typical  multi-model  estimator  is  shown  in 
Fig.  2.  In  this  section,  we  give  a  review  of  the  multi-model  esti¬ 
mation,  with  an  emphasis  on  the  estimate  fusion  strategy. 

Its  first  part  is  composed  of  a  bank  of  parallel  filters  based 
on  different  models.  Each  filter  assimilates  the  data  to  produce 
its  own  estimate.  All  the  estimates  will  then  be  fused  to  give 
the  best  estimate.  Many  options  exist  for  the  elemental  filter, 
such  as  the  KF  for  a  linear  model  or  the  EKF  for  a  nonlinear 
one.  What  is  of  particular  interest  here  is  the  design  of  the 
fusion  strategy. 

Let  us  consider  a  general  system.  Its  unknown  state  at  time 
instant  k  is  denoted  by  xke  [Rnx  and  its  measurement  by  zk<=  Rnz. 
Different  models  are  available  to  describe  the  system,  leading  to 
a  model  set  M  =  {Jt\,  Mi,  •••,  Jt^).  Suppose  that  M\  is  given  by 


Fig.  1.  A  schematic  description  of  adaptive  SoC  estimation. 


Fig.  2.  The  structure  of  a  multi-model  estimator. 
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where  f1  and  h1  are  C1  functions  to  represent  the  state  transition  and 
measurement,  respectively,  and  {wJJ  and  {vJJ  are  uncorrelated, 
zero-mean,  white  Gaussian  noise  sequences  with  covariances  Q!k  > 
0  and  Rlk  >  0,  respectively.  While  assuming  that  the  true  system 
coincides  with  one  model  at  each  time  instant,  we  do  not  know 
which  model  matches  the  system  at  any  time.  Thus  a  probabilistic 
description  is  used.  Let  s/<  denote  the  system  running  status  at  k.  It 
may  take  any  Jt\  for  i  =  1,  2,  ...,  N  to  address  the  uncertainty  of 
model  matching.  The  probability  of  the  event  sk  =  is  denoted  as 
p(sk  =  Ji j),  or  simply,  p{slk).  In  other  words,  p(slk)  indicates  the  a 
priori  probability  that  the  true  model  is  at  time  k.  Obviously, 

EfLiP(4)  =  !• 

From  a  statistical  perspective,  xk  and  zk  are  continuous  random 
variables  and  s/<  a  discrete  one.  Without  causing  confusion,  we  use 
the  symbol  p  to  denote  the  probability  density  function  (pdf), 
probability  mass  function  (pmf)  or  mixed  pdf-pmf  in  the  sequel  for 
convenience.  We  define  the  information  set  as  Zk  =  {z^ ,  z2 , . . . ,  zk} 
and  intend  to  estimate  x/<  from  Zk ,  hence  considering  p(xk \zk).  By 
the  Bayes’  theorem,  we  have 


N 

P(xk|Zk)  =  ^2p(xk,s[  |zk) 

1  =  1 
N 

^p{xk\^zk)p(s[\zky 


(2) 


i  =  l 


When  p(xk\Zk)  becomes  available,  we  can  carry  out  minimum- 
mean-square-error  (MMSE)  estimation  or  Maximum  a  Posteriori 
(MAP)  estimation  of  xk: 


MMSE  : 


:  xk|fc  =  E(xk|Zk)  =  J xklkp(xk\Zk)dxk, 


MAP  :  xk|k  =  argmaxpfxk|Zk). 

Independent  of  the  method  (MMSE  or  MAP)  used,  it  follows 
from  (2)  that 


N 

*k\k  =  £*L|fcp(sk|  Zk), 
1  =  1 


(3) 


where  xk\k  is  the  estimate  of  xk  based  on  the  model  An 
observation  from  this  analysis  is  that  p(s^ \zk)  turns  out  to  be  a 
probabilistic  weight  coefficient.  The  associated  estimation  error 
covariance  is 


■  k\k 


E  [(Xk-Xfe)  (xk-xk)T|zk] 

J  (*k  -  Xk)  (xk  -  Xk)  Tp(xk|Zk)dxk 


p(4|zfe) 


p(4-zfc)  _  p(zki4>zk-i)p(4|zk- 


P(Z  k)  P(Zk|Zk_!) 

p(zki4,zk_i)p(4|zk_1) 
Ef=iP(zk,4|zk_i) 
p(Zki4,zk 1)p(4|zk-1) 


Ef=ip(zki4,zk-i)p(4|zk- 

Furthermore,  we  have 

p(zfcl4-Zk-l)  =  J  P  (zk  >  Xfc  |  Sk ,  Zk_  •[ )  dxk 

P  (zfc|xk,  4-Zk-i)p  (xkl4.  Zk-i )  dxk 
p(zfc  |x,o  4)p  (xk|sk,  zk_i  )dxk. 


(5) 


/ - 
-/< 

Under  the  mildly  simplified  assumption  that 

P(Zfc|xk,s[)  ■^•(h,(xk),R[)andp(xk|s^,zk_1)  =  ^'(xk|k_i,Pk|k_1). 

p(zk\slk,Zk_^)  can  be  approximated  as 


p(zfcl4>zk-i)  ~(2u)  "41  2exp 
where 

4  =  zk  -  h'  (xk|k-l)> 

Sk  =  4pj£|k_1(H')T 

r  (x/<|fc-i)- 


~\[zk 


isir1^ 


■rL, 


Hi 


dW 

ax 


Furthermore, 


\  p(zhI4)p(4)  /a 


since  p(Zfe_1|s|i)  =  1  and  p(Zk_i)  =  1  because  is  an  event 
with  probability  1  at  time  k  -  1.  If  we  define  plk  =  p(slk \Zk)  and 
wlk  =  p(zk\slk,Zk_^)  and  suppose  Tzlk  =  p(slk),  (5)  becomes 


WfTU 


Pk  = 


k“k 


(6) 


Flence,  by  (3)  and  (4),  the  fusion  strategy,  or  the  fuser  as  is  called, 
is  given  by 


xk\k  —  xk\kPk’ 
i  =  l 


i  =  l 
N 

=  I  {xk-  xfc)  (xk  -  Xfc) T p (xk|4,  Zfc)  dxfcp  (si | Zfc) 
1  =  1  J 

N  r 

=  E 

i  =  1  L 


|Zfc)dxk 

N 

r  •  Ti 

1  / 

Pk\k  =  J2 
1=1 

Pk|k  +  (xk  -  Xfc)  (Xfc  -  Xfc) 

Mk- 


(7) 


(8) 


-fc|fc  f  (xfc  -  4)  (xk  -  xi)  P  (4 1  zk)  • 


P//I1/  + 


(4) 

Let  us  take  a  closer  look  at  p(slk\Zk).  Using  the  Bayes’  theorem 
again,  we  see  that 


The  final  conclusion  drawn  from  this  analysis  is  as  follows:  the 
fused  estimate  (covariance)  is  a  linear  weighted  combination  of  the 
estimates  from  the  elemental  filters.  It  can  be  noted  that 

•  The  estimation  is  based  on  a  series  of  elemental  filters  and  the 
fusion.  The  process  is  similar  to  a  ‘  weight-based  reconciliation’, 
which  balances  the  role  that  different  models  play  in  the  esti¬ 
mation  task. 
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•  The  residuals  of  the  elemental  filter  based  on  the  ‘correct’  model 
that  best  match  the  true  system  is  expected  to  be  remarkably 
smaller  than  those  of  the  others  [20].  As  a  result,  the  probabi¬ 
listic  weight  associated  to  this  filter,  say,  will  tend  to  increase 
and  downplay  the  others.  The  fused  estimate  will  approach  the 
estimate  based  on  the  correct  model. 


where 

Xfc  =  [SoCk  I<2  I<3  R] T , 

fi(xi,uk)  =  xi-[KoOOO]Tuk, 

^i(xk,uk)=i<'i+x^2ln(T1  +xkl)+x£3ln(T2  +  l  -x^-x^U);. 
Analogously,  by  fixing  Kj  for  i  =  1,  2,  3,  we  have 


3.  Battery  models  and  observability  analysis 


=  f2(x^ufc), 
=  h2(x2k,uk), 


(13) 


We  investigate  the  battery  modeling  in  this  section.  We  first 
develop  two  submodels  from  a  slightly  modified  Nernst  model 
and  then  analyze  the  local  observability  properties  for  each 
one. 

3.1.  Construction  of  multiple  battery  models 


where 

Xfc  =  [SoCk  I<0  R]r, 
f2(x2,uk):m[x2A-x22uk  x22  x23}T , 

h2(x2Mk)  =  K-l  +/<2ln(T1  +xf0)+K3\n(T2  +  -[-x20)-x2k3uk. 


A  battery  model  consists  of  a  set  of  equations  that  relate  the 
input  (charging/discharging  current),  the  state  variables  (e.g.,  SoC) 
and  the  output  (terminal  voltage).  Various  models  have  been  pro¬ 
posed  and  used,  depending  on  the  specific  purposes.  For  SoC  esti¬ 
mation,  we  consider  the  Nernst  model  here  [5]: 

yk  =  1<\  +  K2ln(SoCk)  +  K3ln(l  -  SoC,<)  -  Rukl  (9) 

where  yk  is  the  terminal  voltage  at  time  instant  k ,  uk  is  the  applied 
current  (u  >  0  for  discharging  and  u  <  0  for  charging),  R  is  the  in¬ 
ternal  resistance,  and  JCZ  for  i  =  1,  2,  3  are  constants.  To  make  (9) 
more  capable  of  grasping  the  dynamics  of  certain  batteries,  we 
propose  the  following  modification: 

yk  =  K]  +  /C2ln(Tj  +  SoC/J  +  K 3ln(i2  +  1  —  SoCk)  —  Ruk ,  (10) 


Remark  1.  In  an  implicit  manner,  Jt\  places  more  confidence  on 
the  state  equation  (11),  assuming  that  K o  is  accurate,  and  the  belief 
in  the  measurement  equation  (10)  is  emphasized  in  Ji2  similarly. 
Nevertheless,  it  is  noteworthy  that  the  confidence  level  on  each 
submodel  during  the  estimation  process  is  dynamically  determined 
by  the  fusion  strategy  outlined  earlier  in  Section  2. 

Remark  2.  An  extended  series  can  be  constructed  on  the  basis  of 
each  submodel  if  we  let  the  parameters  take  different  values  that 
are  believed  to  be  close  or  equal  to  the  truth.  For  instance,  the 
Coulombic  efficiency  may  be  100%,  90%  or  even  80%  depending  on 
the  operating  conditions.  Then  Jt\  will  give  birth  to  three  more 
submodels  if  I<0  assumes  AT/Co,  0.9AT/Q)  and  O.8AT/C0,  respec¬ 
tively.  This  allows  considerable  flexibility  for  us  to  describe  the 
battery  dynamics  and  brings  improvements  to  the  single-model 
case. 


where  two  additional  constants  ti  and  t2  are  added.  In  above, 
K\  + 1<2 ln(ii  +  SoCfc)  +  /C3ln(T2  + 1  -  SoC/<)  in  (10)  can  be  regarded  as 
the  open-circuit  voltage  (OCV)  term.  The  dynamic  change  of  the 
SoC  is  described  by  the  integration  of  the  current  over  time.  In  the 
discrete  time,  it  is  given  by 

c  _  c  _  T 

SoCk  =  SoC0  -  V-7^Ui, 

i  =  0  c° 


3.2.  Observability  analysis 

It  is  well-known  that  state  estimation  requires  a  ‘certain’  kind  of 
observability  of  the  system.  Flence,  we  will  analyze  the  observ¬ 
ability  properties  of  Jt\  and  Jt2  before  proceeding  to  SoC 
estimation. 

Consider  a  general  single-input-single-output  (SISO)  system 


where  p  is  the  Coulombic  efficiency,  Co  the  nominal  capacity  in 
ampere-hour  (Ah),  and  AT  is  the  sampling  period.  An  equivalent 
difference  equation  is 

SoC/<+1  =  SoC  k-K0uk,  (11) 


where  I<o  =  r\-  AT/Co.  We  then  obtain  a  state-space  model  for  bat¬ 
teries  by  putting  together  (10)  and  (11 ).  The  model  state  is  SoC/<  and 
the  parameters  are  I<i  for  i  =  0,...,3  and  R. 

For  adaptive  SoC  estimation,  we  will  perform  simultaneous 
estimation  of  the  SoC  and  the  parameters.  To  obtain  a  locally 
observable  model,  one  or  several  parameters  usually  need  to  be 
fixed  in  order  to  estimate  the  others  and  the  SoC.  A  few  options  may 
exist  regarding  which  parameters  are  fixed.  Based  on  our  experi¬ 
ence  with  the  considered  model,  we  separate  the  parameters  into 
two  sets,  fix  one  set  and  augment  the  state  vector  to  incorporate  the 
SoC  and  the  other  set.  Accordingly,  two  submodels  will  be 
constructed. 

Letting  I<0  and  I<i  be  fixed,  the  first  one  can  be  obtained: 


xLi  =fl(xfc-ufc)> 
yk  =  h\xl,uk), 


(12) 


fXfc+i  =f  (xk,uk), 
\  yk  =  h(xk,uk), 


(14) 


where  xeX  of  dimension  n,  ye  Y  and  ue  U.  We  assume  that  1)  X 
and  Y  connected,  second  countable,  Flausdorff,  differentiable 
manifolds  of  class  C9  with  qe  X,  2)  QJ  is  an  open  interval  of  IR,  and  3) 
f:XxU->X  and  h  :  X  ->  Y  are  of  class  Cf.  For  convenience,  f(x,u)  is 
denoted  as  f'(x),  and  h(f(x,  u0),  iq)  =  hUl  °fUo(x).  Following 
[23,24],  the  local  observability  for  f  is  defined  as  follows: 


Definition  1.  Distinguishability 

Two  states  x  and  x *  are  said  to  be  indistinguishable,  written  as 
x  x*,  if  for  each  l  ^  0  and  for  each  input  sequence,  {u0,  ...,U|}e  V1, 
we  have 


/iu'°fUl-1o...ofU0(x)  =  hUl  ofu'-Y..ofu°(x*). 

Otherwise,  they  are  distinguishable. 

Definition  2.  Local  observability 

The  system  ^  is  locally  observable  if  for  any  state  x°eX,  there 
exists  a  neighborhood  D  of  x°  such  that,x<=>x *  implies  x  =  x*  for  each 

X,X*e  O. 

By  Definitions  1  and  2,  local  observability  means  that  x°  can  be 
distinguished  from  its  neighbors  given  the  input  sequence  {uq,  ..., 
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ui}  and  the  output  sequence  yo,  ...,  y/.  It  should  be  noted  that  this 
definition  of  observability  depends  not  only  on  the  system  itself  but 
also  on  the  applied  inputs,  unlike  the  uniform  observability  for  any 
inputs  defined  in  [25].  While  one  sees  various  definitions  of 
nonlinear  observability  in  the  literature,  this  does  not  obstruct  our 
discussion  since  they  are  usually  about  ‘different  measurements 
results  from  different  initial  states  (for  admissible  inputs)’.  The 
interested  reader  can  refer  to  the  literature  on  the  subject,  e.g.,  [26]. 

To  address  the  observability  condition,  the  following  sets  of 
functions  are  defined: 

=  {&(•)}, 

Qt  =  o o...ofu°(.)  :  Uj-eUVi  =  l,--,jand  1  <j  <  /}, 

C2  Uj>qOj. 


satisfying  the  condition  in  Theorem  1.  Such  a  condition  imposed 
on  the  input  is  a  mild  constraint  that  can  be  easily  satisfied  when 
a  battery  is  in  use. 

•  We  can  analogously  determine  that  Jf2  is  also  locally  observable 
if  a  suitable  {uj J  is  used  to  excite  the  system. 

•  Additional  submodels  other  than  Jt\  and  Jt2  can  be  con¬ 
structed  by  fixing  different  parameters.  An  example  is  to  fix  only 
7<j,  which  will  lead  to  another  locally  observable  model.  How¬ 
ever,  no  matter  how  many  submodels  are  used,  the  essence  of 
multi-model  adaptive  SoC  estimation  remains  the  same,  as  will 
be  seen  in  Section  4.  It  is  also  noteworthy  that  the  resultant 
submodel  will  be  unobservable  if  all  the  parameters  are 
assumed  unknown. 


An  observability  criterion  is  presented  in  the  following  theorem, 
please  see  Ref.  [23]  for  the  proof. 

Theorem  1.  In  Ref.  [23]  If  dim  dQ(x)  =  n  VxeX,  then  the  system 
f  is  locally  observable. 

Theorem  1  gives  a  sufficient  condition  to  determine  the  local 
observability  by  relating  it  to  the  full  dimensionality  of  the  codis¬ 
tribution  dQ.  Now  the  local  observability  of  Jt\  and  Ji2  can  be 
analyzed  using  Theorem  1.  Let  us  take  Jt \  for  an  example  since  the 
analysis  for  both  follows  similar  lines. 

Note  that  f1  and  h1  are  of  class  C00.  Suppose  that  the  initial  state 
is  xj  for  Jt\  and  that  there  are  L  measurements  {yi, . . y^,}.  By  (12), 
xl  is  given  by 

k- 1 

4=xJ-[K o  0  0 

i  =  0 

Hence,  we  have 

hj[(xj )  =  /11 Uk  o  f 11 Ul-'  o  •  •  •  o  f 1; “°  (xj )  =  Kj 

+X0,2ln  ^T1  +  x0,l -Kb&)  +x0,3ln  ^2  +  1  — x0,!  +K0 
—  x0,4U/o 

where  hie  Q.  Define  a  matrix  J  with  dimensions  Lx  4: 


4.  Multi-model  adaptive  SoC  estimation 

In  this  section,  we  study  multi-model  adaptive  SoC  estimation 
on  the  basis  of  Sections  2  and  3.  An  IEKF-based  elemental  filter  will 
be  applied  to  Jt\  and  J(2,  respectively,  for  adaptive  SoC  estimation. 
The  overall  estimate  will  be  obtained  by  fusing  all  the  estimates  for 
the  elemental  filters,  leading  to  the  MM-AdaSoC  algorithm. 

4.1.  Adaptive  SoC  estimation 

Adaptive  SoC  estimation  can  be  attained  via  state  estimation, 
because  the  state  vector  of  each  consists  of  both  the  SoC  variable 
and  the  parameters.  Following  Ref.  [13  ,  we  use  the  IEKF.  As  an 
improved  version  of  the  EKF,  it  is  capable  of  giving  more  accurate 
state  estimates  even  for  highly  nonlinear  systems  by  iteratively 
refining  the  estimate  around  the  current  point  at  each  time  instant. 

Consider  applying  the  IEKF  to  the  system  in  (14).  At  k  -  1,  pre¬ 
diction  can  be  made  about  the  next  time  instant.  The  formulas  are 
as  follows: 

xfc|k-i  =  f(xk-i|k-i.Ufc-i).  (15) 

Pk|k-1  =  Fk-lPfc-l|k-lFk"-i  +Q.>  (15) 

where  x  is  the  estimate  of  x,  P  is  the  error  covariance,  Q,  >  0  is  an 
adjustable  matrix  to  account  for  the  process  noise,  and  F  is  given  by 


[dh] 

d  hi 

ds;iT 

di  \ 

[dxj 

dxj 

dxjJ 

**k- 1  —  ^x/<— 1 1/<— 1 5  Uk-\  J 

The  elements  in  the  k-th  row  of  J  are 


I  _  0 hk  

-  dxf:  ~ 


:0,1  *i+*J,i -Ko  Elio  Ui  *2+l-x;  1+Ko  Elio  Ui 


k- 1 


Jk,2  =i^  =  lnh+xoi-Ko£Ui 

n"  i  =  0 


k- 1 


Jk,3  =  ^  =  ln(  t2  +  1  -x^+Ko^Ui), 


f  _  9h/c  _ 

J/c,4  —  9 —  ~Uk- 

OA0,4 

By  observation,  we  have  the  following  conclusions: 

•  The  submodel  M\  is  locally  observable  if  a  suitable  input 
sequence  {u/J  is  applied.  By  ‘suitable’,  we  mean  that  u/<  varies 
sufficiently  in  magnitude  over  time,  or  in  other  words,  {u/J 
contains  a  rich  mix  of  frequency  contents.  In  this  case,  J  will  have 
full  column  rank,  and  as  a  result,  dim  dQ  has  a  dimension  of  4, 


When  the  measurement  y/<  arrives,  X/<|/<_i  can  be  updated  by  the 
new  information  y^  brings.  The  procedure  is  based  on  iteration.  Let 
£  denote  the  iteration  number  and  xj?k  =  xfe|  fe_1  for  £  =  0.  The 


update  formulas  are 

K®  =  pk|k-1Hj.*-1)  [Hj,£_1)P;{|fe_1Hj,(!_1lT  +R 

yk]  =  h(xkV“1).u/<). 


-1 


-  V 

xk\k  ~  xk\k-l 


|x(®)  L,  _  _  |j(C— 1)  ( C  _  1)\1 

[yi<  yk  Hk  (vx/c|/<l  xk\k  )\ ’ 


where  R  >  0  accounts  for  the  measurement  noise  and 


(«) 


dh 


(17) 

(18) 

(19) 


The  iteration  process  stops  when  £  achieves  the  pre-specified 
maximum  iteration  number  £max  or  when  the  error  between  two 
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Table  1 

The  MM-AdaSoC  algorithm:  Adaptive  SoC  estimation  using  multiple  models. 
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initialize  the  implementation:  k  =  0,  =  Xq,  =  <5T,  where  8l  »  0,  for  /  =  1,2 

repeat 

k^k+l 

IEKF  based  adaptive  SoC  estimation: 


4:  for  i  —  1  to  2  do 

5:  import  the  submodel  Mi 

.Mrbased prediction  (time-update): 

6:  project  the  state  ahead  to  obtain  *luk_x 

7:  project  the  error  covariance  ahead  to  obtain  P^_j 

dt 


Flt-1  -  ^7  (**-11*-!’“*)  »  ^k\k-\  -  FA:-1P£-1|  £-lF£-l  +  Q* 


Mi-based  update  (measurement-update) : 
initialize  the  iteration  procedure:  1  =  0,  x^  =  j 


9 

while  l  <  4iax  do 

10 

£<-£+! 

11 

compute  the  Kalman  gain  matrix 

_  dh‘  fV(<)  u\ 

H*  ~dxi{*k 

_  pi  1) 

K£  -  v^k-xtt-k 

[af 

12 

:  update  the  state  estimate 

Xk\k  Xk\k—\  '  ^k 

\yk-h‘(if'\uk)- 

■Hf 

13 

end  while 

14 

assign  = 

15 

update  the  error  covariance 

P£|£ 

=  [l-K<w>H<w)] 

P£|£- 

16 

:  export  ^ -based  SoC  estimate  SoC^  =  xj^  x 

17; 

end  for 

Estimation  fusion 

18 

determine  the  probability  nl  that  the  battery  runs  on  Mi  for  i  =  1,2  with 

lK  = 

19 

for  i  =  1  to  2  do 

20 

compute  the  initial  weights 

%k-xW  }  +*1 


H1/VriT) 


Hi  =  H  (Vl’“‘)  '  4  =  HiPi|t-l(Hi>T  +R‘ 

>i|jt-i  =^-^-1,  4  =  (27T)-t(4)-2exp 


Ml! 

25i 


21:  end  for 

22:  compute  the  normalized  weights 


23:  fuse  the  SoC  estimates  from  M\  and  M2 


24:  until  SoC  estimation  task  ends 


Pk  = 


_  2  _ . 

SoC*  =  X  SoC£Mi 


for  i  =  1,2 


consecutive  iterations  is  less  than  the  pre-selected  tolerance  level. 
Then  xk\k  =  xj^H,  and  the  associated  error  covariance  is  given  by 

pk\k  =  [i-<max)<max)]p,<i,<-i- 

Following  the  above  description,  the  IEKF  can  be  applied  as  an 
elemental  filter  to  Jt\  and  Jt 2.  The  resultant  state  estimates  are 

xl\k  and  xj^,  respectively.  Accordingly,  the  SoC  estimates  are 

—  i  — - — 2  _2 

denoted  as  SoCk  =  l  and  SoCk  =  xk|k  l,  respectively. 

4.2.  MM-AdaSoc:  multi-model  adaptive  SoC  estimation 

- - -1  - - -2 

The  SoC  estimates  produced  from  M\  and  SoC,<  and  SoCfc, 
respectively,  can  be  combined  weightedly  to  generate  the  overall 
estimate  SoCk.  In  the  light  of  the  fusion  strategy  in  (7)  and  (8),  we 
have 


sdck  =  f)  2S SCfcMi,  (20) 

i  =  1 

where  the  weight  coefficient  m  for  i  =  1, 2  can  be  determined  using 

(6). 

Putting  together  the  results,  we  obtain  the  MM-AdaSoc  algo¬ 
rithm,  as  is  summarized  in  Table  1. 

Remark  3.  The  underlying  idea  of  the  proposed  MM-AdaSoc  al¬ 
gorithm  is  that  the  IEKF-based  adaptive  SoC  estimation  is  carried  out 
for  multiple  models  and  then  the  estimation  results  are  fused  to  yield 
the  overall  SoC  estimate.  For  the  MM-AdaSoc,  the  recursive  and  real¬ 
time  implementation  cuts  down  the  amount  of  stored  data.  Mean¬ 
while,  higher  estimation  accuracy  is  achieved,  because  the  update 
procedure  relies  on  iterative  searching  at  each  recursion.  Another 
noteworthy  advantage  is  that  a  good  balance  is  maintained  between 
the  estimation  performance  and  the  computational  complexity, 
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conceding  a  generally  linear  moderate  increase  of  the  demanded 
computing  power  depending  on  the  number  of  models  used. 

Remark  4.  The  applicability  of  the  proposed  MM-AdaSoC  algo¬ 
rithm  to  different  types  of  batteries  is  quite  promising.  Due  to  its 
parameterized  characterization,  the  Nernst  model  has  been  found 
capable  of  describing  the  dynamics  of  many  batteries,  e.g.,  nickel 
metal  hydride  (NiMH),  LiMn204  and  UC0O2.  As  a  result,  the  mm- 
AdaSoC  algorithm  can  be  well  applied  to  such  batteries  for  its 
construction  based  on  the  Nernst  model. 

Remark  5.  Not  limited  to  the  MM-AdaSoC  algorithm  at  all,  the  role 
that  multi-model  estimation  can  play  is  more  profound.  It  can  be 
developed  as  a  framework,  within  which  a  variety  of  advanced 
estimation  methods  can  be  built  for  battery  applications.  Here,  we 
identify  five  potential  sources  of  multiple  models: 


Here,  ko  =  1.0683  x  107  per  hour,  F^s/Rg  =  6.789  K,  the  current 
efficiency  7]  =  1  for  discharging  and  0.99  for  charging  near  50%  SoC, 
the  nominal  capacity  Co  =  1.25  Ah,  the  reference  temperature 
T  =  35  °C  (308.15  K),  and  the  sampling  period  AT  =  1  s.  The  initial 
SoC  is  assumed  to  be  50%.  The  terminal  voltage  is  equal  to  the  OCV 
plus  internal-resistance-induced  drop,  that  is, 


y/<  —  ^oc,/<  R^k- 

The  OCV  Voc  is  given  by  the  following  equation  with  the  inclu¬ 
sion  of  voltage  hysteresis: 


Voc  =  Un 


neF 


SoCk  -  n 
1  -  SoQ 


+  V; 


H.lo 


where  the  varying  voltage  hysteresis  VH  is  characterized  by  an 
empirical  expression 


•  a  set  of  submodels  established  from  a  battery  model  by  fixing 
certain  parameters  for  adaptive  SoC  estimation,  as  we  have  done 
in  this  paper, 

•  a  set  of  submodels  established  from  a  model  by  assuming 
different  sets  of  values  for  model  parameters, 

•  a  set  of  different  models  constructed  in  different  ways,  such  as 
an  equivalent-circuit  model  and  an  electrochemical-principles- 
based  model, 

•  a  set  of  models  capturing  different  characteristics  of  batteries, 
e.g.,  the  charging  and  discharging  processes,  cycling  and  aging 
effects,  and 

•  a  multitude  of  (sub)models  combining  the  above  four  cases. 
The  multi-model  approach  promises  three-fold  benefits. 

•  It  better  apprehends  the  battery  dynamics  known  to  be  complex 
and  multi-faceted,  thus  promoting  the  accuracy  and  robustness 
of  SoC  estimation. 

•  It  reduces  the  complexity  of  estimator  design,  especially  when 
highly  nonlinear  battery  dynamics  are  involved,  in  a  ‘divide- 
and-conquer’  manner.  Simple  and  elegant  solutions  will  be 
achieved  and  theoretical  analysis  is  made  easier. 

•  It  can  even  provide  useful  model  interpretation  and  comparison 
in  some  circumstances. 

To  fully  realize  its  potential  and  benefits,  multi-model  estimation/ 
control  for  batteries  needs  to  be  further  studied  in  the  future. 


5.  Application  examples 


In  this  section,  we  present  two  examples  using  simulation  and 
experiment  data,  respectively,  to  evaluate  the  MM-AdaSoC 
algorithm. 


Example  1.  This  example  is  based  on  simulation  with  a  model 
used  in  [15]  for  a  NiMH  battery  system.  For  simulation  purpose,  we 
employ  certain  minor  modification,  but  the  obtained  model  is  still  a 
sufficiently  accurate  representation  of  the  NiMH  battery  dynamics 
in  most  circumstances.  The  change  of  SoC  is  governed  by 


S°Ck+1  =  SoQ  -  ■ 


■  AT  sD(jref)-AT 


-uk- 


100 


where  the  third  term  on  the  right-hand  side  represents  self¬ 
discharge  with 


So(T'ref)  =  fcoexp 


Eas  ) 

^g^ref  J 


SoC. 


Vr/<+i  -  VH)/<  -P-rj- [Vu,max  +  sign(U/J-VH/J -AF-i^. 

The  resistance  R  is  described  by 

n 

R  =  ^QjSod. 
j  =  0 

The  Farady’s  constant  F  =  9.6485 xlCT4  C  mol-1,  ne  =  1, 
Rg  =  8.314  J  mol-1  K-1,  U0  =  1.37  V,  n  =  0.08,  0  =  3  x  10-5  C-1, 

Vh ,max  =  0.05  V.  The  initial  Vh  is  0.005  V.  In  addition, 

a0  '=  4.1252  x  10-2,  ax  =  8.9691  x  10-4,  a2  =  1.6760  x  10-5, 

a3  =  -1.4435  x  10-7  and  a4  =  4.7223  x  10  10. 

During  the  simulation,  we  do  not  assume  that  this  model  is  fully 
available  for  SoC  estimation.  Instead,  the  modified  Nernst  model 
presented  in  Section  3  will  be  used  for  approximate  description  of 
the  above  true  model.  Let  K 0  =  r/  •  AT/Co,  K\  =  Do,  K 2  =  I<3  =  RgTref/ 
(neF),  ii  =  -n,  t2  =  0  and  R  =  do.  The  current  signal  applied  as  the 
input  to  the  battery  was  a  pseudo-random  binary  sequence  (PRBS) 
stretched  by  100  times  over  the  time  axis.  Its  magnitude  is  1  A.  A 
view  of  the  input  current  and  output  voltage  during  the  first  2000 
time  instants  is  given  in  Fig.  3(a).  Let  the  true  initial  SoC  be  55%  and 
the  initial  SoC  estimate  be  65%.  The  initial  weights  assigned  to  M\ 
and  Jt2  are  0.7  and  0.3,  respectively. 

In  this  setting,  we  face  hysteresis,  model  mismatch  and  incorrect 
initial  estimate,  which  together  make  SoC  estimation  a  tougher 
challenge.  As  described  in  previous  sections,  we  consider  two 
submodels,  with  the  first  one  assuming  known  K 0  and  7<j  and  the 
second  assuming  known  Kj,  I<2  and  K3.  When  the  MM-AdaSoC  al¬ 
gorithm  is  applied,  Fig.  3(b)  shows  the  estimation  of  the  SoC  over 
time.  We  see  that  both  Jt\  -  and  Jt2- based  estimates  differ  from  the 
actual  values  with  bounded  errors.  However,  it  turns  out  that  the 
overall  estimates  given  by  the  MM-AdaSoC  algorithm  become  more 
accurate,  demonstrating  that  the  estimation  errors  can  be  reduced 
effectively  by  fusion  of  the  multi-model  estimates.  The  weights  of 
the  two  models  are  compared  in  Fig.  3(c).  Obviously,  Jt2  weighs 
much  more  than  M\  in  this  case.  This  is  because  M\  is  sensitive  to 
the  initial  SoC  estimate,  relying  on  the  state  equation  based  on 
Coulomb  counting.  Furthermore,  we  also  apply  the  well-known 
EKF  to  the  modified  Nernst  model  with  known  I<i  for  i  =  0, 1,  2,  3 
for  SoC  estimation.  As  shown  in  Fig.  3(d),  the  EKF  yields  unreliable 
results  in  this  situation  in  comparison  to  the  MM-AdaSoC  algorithm. 

Example  2.  For  the  experimental  evaluation  of  the  MM-AdaSoC 
algorithm,  data  was  collected  from  a  LiMn204/hard-carbon  battery 
in  the  Advanced  Technology  R&D  Center,  Mitsubishi  Electric  Cor¬ 
poration.  The  experiment  was  conducted  using  a  rechargeable 
battery  test  equipment  produced  by  Fujitsu  Telecom  Networks.  The 
current  input  was  a  PRBS  signal  stretched  by  10  times  over  the  time 
axis  with  a  magnitude  of  5  A.  Despite  many  other  options,  we  chose 
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(a) 


(b) 


(c) 


(d) 


Fig.  3.  (a)  The  input  current-output  voltage  profile;  (b)  SoC  estimates  versus  truth  over  time;  (c)  fusion  weights  for  M\  and  M2  versus  time;  (d)  comparison  with  the  EKF. 


the  PRBS  because  it  has  white-noise-like  properties  and  is  admis¬ 
sible  for  observability.  The  profile  of  the  input  current  and  the 
output  voltage  is  shown  in  Fig.  4.  The  battery  has  a  nominal  ca¬ 
pacity  of  4.93  Ah.  The  sampling  period  was  1  s.  During  the  exper¬ 
iment,  the  ambient  temperature  in  the  chamber  was  maintained  at 
25.8  °C. 


We  consider  the  model  in  (10)  and  (11).  The  Coulombic  effi¬ 
ciency  constant  I<0  =  5.6342  x  10  5  when  u/<  >  0  (100%  for  dis¬ 
charging)  and  K0  =  4.7891  x  10~5  when  u/<  <  0  (85%  for  charging). 
From  the  SoC-OCV  data  collected  from  this  type  of  batteries,  it  can 
be  determined  that  K \  =  1.294,  I<2  =  0.0984,  I<3  =  3.972, 
Tl  =  I2  =  0.3. 


Fig.  4.  The  input-output  (current-voltage)  profile. 
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Fig.  5.  (a)  SoC  estimates  versus  time;  (b)  SoC  estimates  during  the  initial  450  s;  (c)  fusion  weights  for  Jt\  and  Jt2  versus  time;  (d)  fusion  weights  during  the  initial  450  s;  (e) 
comparison  between  the  true  and  the  one-step-forward  predicted  voltage. 
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As  aforementioned,  the  actual  values  of  the  parameters  I<i  for 
i  =  0, . . 3  can  change  as  a  result  of  the  operating  conditions.  Hence, 
rather  than  depending  fully  on  their  nominal  values,  we  perform 
multi-model  adaptive  SoC  estimation  by  applying  the  MM-AdaSoC 
algorithm.  The  construction  of  two  submodels  from  (10)  and  (11 )  is 
described  in  Section  3.1. 

The  SoC  estimation  results  are  shown  in  Fig.  5.  The  full  view  over 
the  available  experimental  data  is  given  in  Fig.  5(a).  The  initial  SoC 
of  the  battery  is  known  to  be  approximately  50%.  It  is  seen  that 
there  is  a  difference  of  approximately  5%  between  the  Ji\  -based 
and  Jt2 -based  estimates.  Based  on  our  experience,  M\  tends  to 
yield  conservative  estimates  in  this  case  and  Jt 2  does  the  opposite. 
The  MM-AdaSoC  algorithm,  through  the  fusion  strategy,  makes 
adjustment  to  give  neutralized  overall  estimates.  Although  the  true 
SoC  data  are  not  available,  we  still  judge  that  the  estimates  are  close 
to  the  truth,  based  on  our  a  priori  knowledge  about  the  battery 
behavior.  Fig.  5(b)  illustrates  what  happens  during  the  initial  450  s. 
It  is  seen  from  Fig.  5(a)  and  (b)  that  the  overall  estimates  are  closer 
to  those  based  on  Ji2 .  This  is  verified  in  Fig.  5(c)  and  (d),  where  the 
weight  for  Jt\  fluctuates  slightly  around  0.63  and  p2  around  0.37. 
Thus,  with  a  larger  weight,  Jt\  is  given  more  confidence  than  Jl2 
by  the  MM-AdaSoc  algorithm  during  the  implementation.  It  is  un¬ 
derstood  that  the  fusion  depends  on  the  performance  of  one-step- 
forward  prediction  of  the  terminal  voltage.  Fig.  5(e)  compares  the 
measured  data  with  the  prediction  based  on  Jt\  and  Jt2,  respec¬ 
tively.  The  prediction  is  satisfactory  for  both  submodels,  but  Ji2  is 
observed  to  lead  to  the  better  predicted  voltage. 

From  the  above  results,  we  believe  that  the  MM-AdaSoc  algo¬ 
rithm  is  quite  effective,  supported  by  the  findings  that  the  obtained 
SoC  estimates  exhibit  considerable  accuracy  and  that  the  voltage 
prediction  approximates  the  truth  well.  Through  experiments  with 
charging/discharging  rates  of  0.5  A,  1  A,  10  A  and  15  A,  we  consis¬ 
tently  observe  similar  estimation  results,  which  shows  that  the 
applicability  of  the  model  and  the  power  of  the  MM-AdaSoc 
algorithm. 

6.  Conclusions 

Development  of  adaptive  approaches  for  SoC  estimation  is  of 
practical  significance,  because  battery  dynamics  are  often  hard  to 
fully  determine  and  are  time-varying.  We  are  focused  on  improving 
the  adaptive  SoC  estimation  via  launching  a  multi-model  strategy 
in  this  paper,  motivated  by  the  proven  success  of  multi-model 
estimation  in  addressing  problems  involving  structural  and 
parameter  changes. 

The  main  contribution  of  this  paper  is  the  construction  and 
validation  of  the  MM-AdaSoc  algorithm.  It  is  built  to  estimate  a 
battery’s  SoC  in  real  time  through  carrying  out  simultaneous  state 
and  parameter  estimation  on  a  set  of  (sub)models.  We  first 
construct  two  submodels  from  a  general  state-space  battery  model 
by  fixing  different  parameters,  with  both  shown  to  be  locally 


observable  with  admissible  inputs.  The  well-known  IEKF  is  then 
applied  to  each  submodel  to  produce  the  SoC  and  parameter  esti¬ 
mates.  The  final  overall  estimates  are  generated  by  fusing  the 
submodel-based  estimates,  and  it  is  shown  that  the  fusion  is  a 
linear  weighted  combination  of  the  estimates.  Simulation  and 
experimental  results  are  presented  to  demonstrate  and  validate  the 
effectiveness  of  the  algorithm. 

Apart  from  the  MM-AdaSoc  algorithm,  we  also  emphasize  the 
potential  of  the  multi-model  framework  for  battery  applications. 
The  initial  success  reported  in  this  paper  would  provide  strong 
incentives  for  further  development  of  a  wide  range  of  methods 
based  on  multiple  models  to  better  monitor  the  status  and  health  of 
a  battery. 
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